Quantum phase measurement and Gauss sum factorization of large integers in a 
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We study the implementation of quantum phase measurement in a superconducting circuit, where 
two Josephson phase qubits are coupled to the photon field inside a resonator. We show that the 
relative phase of the superposition of two Fock states can be imprinted in one of the qubits. The 
qubit can thus be used to probe and store the quantum coherence of two distinguishable Fock states 
of the single-mode photon field inside the resonator. The effects of dissipation of the photon field 
on the phase detection are investigated. We find that the visibilities can be greatly enhanced if the 
Kerr nonlinearity is exploited. We also show that the phase measurement method can be used to 
perform the Gauss sum factorization of numbers (> 10 4 ) into a product of prime integers, as well 
as to precisely measure both the resonator's frequency and the nonlinear interaction strength. The 
largest factorizable number is mainly limited by the coherence time. If the relaxation time of the 
resonator were to be ~ 10 fis (~ 1 ms), then the largest factorizable number can be > 10 4 iV (> 10 7 7V), 
where N is the number of photons in the resonator. 

PACS numbers: 03.65. Vf, 42.50.Dv, 85.25.Cp 
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I. INTRODUCTION 

The superposition of states is a fundamental feature 
of quantum mechanics. Recently, the arbitrary superpo- 
sition of Fock states [l|, [1] has been produced in a su- 
perconducting resonator with a Josephson phase qubit 
[3| . This offers novel ways to directly study the quantum 
coherence of the photon field, i.e., superposition of num- 
ber states. Also, this strongly coupled qubit-resonator 
system 0, [E], may be useful for quantum information 
processing (QIP). 



A. Quantum phase measurement 

We theoretically study the quantum phase measure- 
ment of the photon field in a superconducting resonator 
coupled to two phase qubits. We consider probing the 
quantum coherence of the superposition state 



\* N ) = ^[\0) 



■exp(i<p N )\N)]. 



(1) 



where |0) and \N) are the vacuum and the multi-photon 
state, respectively, and tpjj is the relative phase. This su- 
perposition of states leads to interference fringes. How- 
ever, in a dissipative environment, the quantum coher- 
ence of the superposition in Eq. ([T]) decays rapidly as 
N grows large @. Systematically studying such super- 
positions should provide a better understanding of the 
decoherence process @, 0] ■ 

To measure the quantum state of the photon field, 
Wigner tomography can be used The relative phase 

between two number states can only cause a rotation in 
Wigner phase space without changing the shape of the 
Wigner function Q. Alternatively, here we propose a 
method to transfer the phase information of the photon 



field to the qubit, such that we can determine the relative 
phase precisely by measuring the quantum state of the 
qubit. In this way, the qubit can be used to store and de- 
tect the quantum coherence of an arbitrary superposition 
of two photon number states. 

The superposition of the vacuum and the single-photon 
state l^i) for N = 1 in Eq. (Q]) in a microwave cavity has 
been used for the quantum memory of an atomic qubit 
0. This experiment has demonstrated that the quan- 
tum information of the qubit can be transferred to the 
photon field. However, here we find that it is necessary 
to use two qubits to transfer the phase information of the 
superposition of the vacuum and the multi-photon state 
l^jv) i n Eq. ([TJ. The first qubit is used for storing the 
quantum information of the photon field, whereas the 
second qubit is used as an auxiliary qubit to disentan- 
gle the first qubit from the resonator [lOj by repeatedly 
appl ying a controlled-NOT (CNOT) quantum gate (see 
e.g., [iMl). Several proposals (e.g., GSII3) 

have been 

made to implement CNOT gates using superconducting 
qubits. 

Our proposed quantum phase detection method can 
measure the degree of quantum coherence of the photon 
field, which can be determined by the visibility of the 
detection signal [6;]. However, the visibilities are greatly 
reduced due to decoherence, and depend on the quality 
factor [l7j | (i.e., ratio of the frequency and the damping 
rate of the resonator). We find that the visibilities can 
be greatly enhanced if the Kerr nonlinearity of the cavity 
mode (e.g., (l8l.[l9j) is exploited. Thus, the interference of 
the superposition of multi-photon states can be observed 
clearly even in the presence of a dissipative environment. 
Notably, the production of extremely strong Kerr non- 
linear strength via coupled Cooper Pair Boxes (CPBs) 
(20l | using the effect of electromagnetically induced trans- 
parency (EIT) (e.g., [2H - I23I I) has recently been proposed 
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(24 |25| . The EIT effect in a lossiess medium, such as 
quantum dots embedded in a soiid-state substrate, has 
recently been studied |26l |. 



B. Gauss sums 

The big problem with Shor's algorithm [13, [28| is that 
it is very difficult to implement physically for numbers 
larger than 15. Thus, here we use the Gauss sum ap- 
proach [29U331 ] because this is implementable and could 
provide a very valuable test-bed or stepping stone for 
more powerful future implementations. Our proposed 
phase detection scheme can be used for implementing the 
Gauss sum, which can find the factors of a number using 
the periodicity of the sum of the quadratic phase factors 
[29l | . The Gauss sum has been realized with NMR [30j . 
cold atoms [U H3 and using short laser pulses [33| . Here 
we study the implementation of the Gauss sum in super- 
conducting circuits. Using our proposed superconducting 
circuits, factors of integers ~ 10 3 should be obtainable 
by the Gauss sum if the relaxation times were to be sev- 
eral fis. The size of the factors are mainly limited by the 
coherence time of the photon field in the superconduct- 
ing resonator. Thus, the factors of much larger numbers 
could be obtained in the future, when coherence times 
improve. 



C. Measuring the frequency 

We also show that the quantum phase measurement 
approach proposed here can be applied to precision mea- 
surements. This approach enables to precisely determine 
the frequency of the resonator and the strength of the 
Kerr nonlinearity. We show that the superposition of 
multi-photon states can increase the accuracy of precision 
measurements. This can act as a "frequency standard" 
for the other qubits in the circuit. 

This paper is organized as follows: in Sec. II, we de- 
scribe the model of the system studied. In Sec. Ill, we 
present a method to detect the relative phase of the su- 
perposition of the vacuum and the multi-photon state in 
Eq. ([1]). In Sec. IV, we show that the phase detection 
method can be applied to both the Gauss sum for factor- 
ization and precision measurements. We close the paper 
with a summary. In appendix A, we study the interac- 
tion between a qubit and a resonator in the far-detuning 
regime. In appendix B, we discuss the effect of imperfect 
CNOT gate operations on the disentanglement process. 



II. SYSTEM 

We consider two Josephson phase qubits capacitively 
coupled to a superconducting resonator [34| as shown in 
Fig. [T] The Hamiltonian of the qubit-resonator system 



can be written as Q (H = 1) 

2 

H = H Iea + ^2 [-^qbit + ^qbit-res + ^drive] ' 
2 

• [exp (—i<f>j — iuj q jt)aj + + H.c.] j, 



(2) 



uia)a - 



a,- 



(3) 



where a is the annihilation operator of the photon field, 
<jj± and <jj z are the transition and population (Pauli z) 
operators of the qubit j respectively (j = 1,2). Here 
tooj and to are the frequencies of the qubit j and the res- 
onator, respectively. The parameters gj act as coupling 
strengths between the photon field and the qubits. Also, 
Clj and 4>j are the amplitude and the phase of the cou- 
pling between the ground state \g) and the excited state 
|e) in the j-th qubit. The frequency to q j is the frequency 
of the microwave drive of the qubit j. The qubit j can 
be accurately controlled by adjusting the frequency u>oj 
and the parameter Clj via a classical signal 0, 0] • 
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FIG. 1. (Color online) Circuit diagram of the qubits and the 
resonator for quantum phase measurements. The Josephson 
qubits are capacitively coupled to the coplanar waveguide res- 
onator with couplings gi and t/2- The frequency of the qubit 
can be adjusted by the flux bias pulse and can be measured 
by a superconducting quantum interference device (SQUID) 
[34|] . Here f2i and Q2 are the amplitudes of the couplings 
between the ground and excited states. 

It is convenient to work in the interaction picture. By 
considering the unitary transformation 



ttu\ r f t , w oi . w 2 \ 

U [i) = exp ^ - i[uja [ a + — G iz + —^2zji 



the transformed Hamiltonians are 



H^ = UHt)H^ ve U(t), 



(4) 



(5) 



iAjt)cr j+ + H.c] , (6) 
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and 



where 



and 



H. 



qbit— res 

U(t), 

gj [exp(iAjt)a(Tj + + H.c] , 



(7) 

(8) 

(9) 
(10) 



Aj = LOQj — LJ. 

Hereafter, the two transformed Hamiltonians and 
H2 at resonance (i.e., when Aj = and Aj = 0), 

H U) = _^ [ exp (- i(f)j )a j+ + H.c.] , (11) 

(12) 



H 2 ] = 9j(aa j+ +CTj_a t ), 



will be used frequently. Also, the time evolution opera- 
tors U¥\t) = exp[-iH[ j h] and £/ 2 (i) (t) = exp[-iH^t] 
can be obtained explicitly as [23| 



U[ J) (t) = cos f^j 1 + i sin f^j [ cxp (-i^)a j+ 
+H.c], (13) 

and 

U^it) = cos (gjtVaa^)\e)(e\ + cos (gjtVaJa) \g)(g\ 



sin (gdy aJa) 
-la '-^== -a 



3+ 




(14) 



where 1 is the unit operator. 

To transfer the phase information of the photon field to 
the qubit, it is required to switch-on and off the interac- 
tion between the qubit and the resonator. The coupling 
between the qubit and the resonator is fixed, but the 
frequency of the qubit can be adjusted by the bias cur- 
rent [1, H[ • The qubit-resonator interaction can thus be 
turned-off by far-tuning the frequency woj of the phase 
qubit [1, 3- I n appendix A, we give a detailed discus- 
sion of the qubit-resonator coupling in the far-detuning 
regime. 



III. QUANTUM PHASE MEASUREMENT 

We now present a procedure to use the qubits to probe 
and store the quantum coherence of the superposition of 
two Fock states. The relative phase between the super- 
position of the vacuum and the smg/e-photon state in 
Eq. (JT]) can be completely transferred to the qubit such 
that the qubit can act as a probe of the quantum coher- 
ence of the photon field @ . The phase can be determined 
by measuring the qubit's state. 

We require two qubits to measure the phase of the 
superposition state of the vacuum and the mw^i-photon 



state in Eq. (fTJ . One qubit is used for probing the phase 
and the other qubit is used to disentangle the qubit from 
the resonator. In the following subsections, we will de- 
scribe the different schemes for the single- and multi- 
photon cases. 



A. Single-photon case: jO) + exp(iy>i)|l) 

Let us now consider the quantum phase measurement 
of the superposition, [|0) + exp(it/?i)|l)]/v / 2, of the vac- 
uum and the single-photon Fock state in the resonator. 
We can create the superposition state \^>i) in Eq. (JTJ) by 
just using one qubit (2. •''>. 9]. Initially, the product state 
of the vacuum of the resonator and the ground state of 
the qubit 1 is prepared, i.e., | qubit, resonator) = |<?,0). 
We first produce an equal superposition of the states \g) 
and |e) of the qubit 1 by applying a 7r/2-pulse to the 
qubit 1 (i.e., turning on the drive of the qubit 1 for a 
time T — 7r/257i). By applying the time-evolution oper- 
ators U{t)u[ 1] (t) in Eqs. (0} and (H3J) to the state \g), 
we have 



= tf (*)[«»( — ) \g) + 1 exp(-?</>i) sin 

Now consider u[ l \t = T = 7r/2Qi), then 
1 

f i exp(-i4>i)\e}] , 



(15) 



7T xp ( 



iujQit 



*exp(-i#)|e)], (16) 



where 



1*1(0)) 



+ ojQit. The states can thus be written as 
(WP(t = 77/20015,0), (17) 
= -L cxp (^) [|s) + i exp (-i#)|e>] |0>. 

(18) 

Now for U(t = 7r/2fti), then 
1 



l*i(0)) 



fiuj 01 n 



[]<?)+; exp (-#i)|e)]|0). 

(19) 



Next, we turn on the qubit-resonator interaction for a 
time t*, 



7T 
25? 



(20) 



The energy of excited state of qubit 1 then lowers one 
level, to its ground state, and a photon is created in the 
resonator. This can be derived by applying the evolution 
operator U(tl)U^\t$) in Eqs. f|4TT4]t giving 

l^)[/ 2 (1) (^)|e,0} = -iexp(-iutl)\g,l), (21) 
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and the ground state \g)\0) has not changed, 

U(tt)U^\tt)\g,0} = \9,0}. (22) 

Combining Eqs. (|18l21|22j) . we obtain the state 

|*i(ii)) = ^(tl)t^ 1) (tl)|*i(0)) 1 (23) 
= i= | 5 ) { |0) + exp + art*)] 1 1) } . (24) 

Here we have ignored the global phase factor 
exp(iWoi7r/4fii) in Eq. ([19]). 

We then switch-off the qubit-resonator interaction and 
let the system evolve freely for a short period r, such 
that a relative phase is acquired between the two number 
states |0) and |1). The total state, at the time t' = t\ + r, 
becomes 

|$i(0) =eotp(-iwoW)|$i(iI)), (25) 
= -L \g) { |0) + exp -iw(ij +r)]|l)}. 

(26) 

We can now transfer the phase information to the qubit 
1 by switching-on the qubit-resonator interaction for the 
time i* in Eq. (gDJ. The state will swap to |e)|0), 

while the ground state |<?)|0) remains unchanged. The 
state now reads 

|*!(t)> = i={| 9 >-iexpH(^ + ^)]|e>}|0), (27) 



1. Dissipation in the photon field 

The photon field inevitably suffers from the dissipation 
present in realistic situations. The thermal average pho- 
ton number is about zero (~ 1CP 6 ) for a high-frequency 
resonator (~ 40 GHz) at low temperatures (~ 25 mK) 
3]. The time evolution of the density matrix p of the 
photon field can be described by the master equation 

p = —iuj[a)a, p] + T(2apa) — a'ap — pa 1 a), (30) 

where T is the damping rate. We assume that the number 
of thermal photons is negligible. Here we ignore the de- 
coherence effect during the qubit and qubit-resonator op- 
erations because their time durations (T and t* ~ ns in 
3) are extremely short compared to the dissipation time 
r _1 (several ps in Q). Therefore, we consider the dissi- 
pation of the photon field during the free time-evolution 
t. The density matrix p of the photon field at the time 
t can then be found as [35[ 

p=i{2|0)(0|-exp(-2rY)(|0)(0|-|l)(l|) 

+ exp(-rr)[exp(^T + i0 , 1 )|O)(l| +H.c.]}. (31) 

We then follow the same procedures discussed above. 
The probability P e that the qubit is in its excited state 
can be readily obtained 

P e = i[l-e X p(-rr)cos^ 1 ], (32) 
where tp± = lot. 



where t = 2t\ + r and ip± = tot ~ lot. Here we have 



assumed that the period r is much greater than the time 
t*. The relative phase information between the two num- 
ber states is now imprinted in the qubit 1 and the qubit 
1 is disentangled from the photon field in the supercon- 
ducting resonator. 

Let us apply a 7r/2-pulse to the qubit 1 so that the 
final state can be written as 



I* 



1 (t f )) « i e xp(^){[l + e Xp (-V 1 )]| 5 ) 

+^ exp(-i#) [1 - cxpH^)] |e) } |0) , (28) 



where tf = t + T, ip\ = ip.\ — loqiT/2 and T — 7r/2f2i. 
We have also assumed that the time T — 7r/2f2i is much 
shorter than the time t and we approximate ip' ± ~ ip%. 
The excited state of the qubit 1 (with a much higher 
tunneling rate than that of the ground state) can now be 
measured. This can be done by applying a measurement 
pulse and read-out by a SQUID 34]. The phase factor 
can thus be determined from the probability of the ex- 
cited state, which is 



P e = ~(1 - cosipi). 



(29) 



as 



2. Visibility 

The visibility of the quantum coherence can be defined 

Cmax Cmin 



V = 



(33) 



where C max and C m - m denote the maximum and mini- 
mum values of the coherence factor C(r), 



C(t) = exp (—IV) cosipi, 



(34) 



which characterizes the quantum coherence of the super- 
position state. Note that the visibility V is unity when 
dissipation is absent. The visibility V can be obtained as 



V 



1 



r 2 



exp (-rT TO ) 



where 



— arccos 
us 



r 2 



(35) 



(36) 



and 7r/2 < u>T m < 7T. Clearly, a higher visibility can be 
obtained for larger resonator quality factors (i.e., ratios 
of lu and T) . 
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B. Multi-photon case: 0) + exp(i<pu)\N) 
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FIG. 2. (Color online) Circuit diagram of the coupled qubits 
and the resonator. The Josephson qubits can be coupled to 
each other via a capacitor (shown in dark red line). Here A is 
the coupling strength between the two qubits. 

We now investigate the quantum phase detection of the 
superposition of the vacuum and multi-photon states in 
Eq. where N > 2. This superposition state |^jv) 
can be generated by applying an appropriate sequence 
[l| of evolution operators 



U£'(T*)U^(T j )...U^(Tt)U^(n) 



(i) 



(37) 



to the initial state \g, 0), where n and r* are the i-th time 

steps for the time evolution operators and 17%, re- 
spectively, and j is the total number of steps. We assume 
that we can adjust the time duration of the interactions 
in each step. The time steps 73 and t* can be found 
inversely [l| by applying 



U- 



(i)t 



.a 



(38) 



to the state \ty N ) in Eq. JIJ. 

We now consider two qubits (qubit 1 and qubit 2) for 
the quantum measurement. We assume that both qubits 
have the same form, Eq. ©, of the interaction with the 
resonator. Qubit 1 is used as an auxiliary qubit to disen- 
tangle qubit 2 from the resonator. A CNOT gate U c can 
be applied to the two qubits. Details for implementing a 
CNOT gate can be found elsewhere (e.g., Refs. [l5l[T6j). 

Note that we now need to add a qubit-qubit coupling 
term in the Hamiltonian in Eq. (see red capacitor in 
Fig. [2|). The qubits in Fig. Q] can be coupled [M HI S3 
to each other via a capacitor so that "spin-spin" interac- 
tions can be produced. This qubit-qubit coupling Hamil- 
tonian, with the coupling strength A, can be written as 



-ffqb-qb = A((7l + (72- + Ci_<7 2 -| 



(39) 



Alternatively, the qubits can also be coupled via a high- 
excitation-energy quantum circuit such as a supercon- 



ducting resonator [38 40]. A CNOT gate can be real- 
ized, e.g., by either applying an external field 15], pulses 
(l6j or applying iSWAP gates [turning-on the qubit-qubit 
interaction in Eq. (|3"§|) for the time t = 7r/4A] appropri- 
ately [!l|]. By adjusting the frequencies of the two qubits, 
the CNOT gate and the qubit-resonator interaction can 
be independently operated. We choose qubit 1 as the 
control qubit and qubit 2 as the target qubit such that: 
I.9SMSS); \ge)->\ge)', M-Hee); |ee)-s-|eg). 

Let us first prepare the resonator and the two qubits 
in a product state as 



I qubit 1, qubit2, resonator) : 



*iv(0)), 
|0> 



\gg) 



I AO 



V2 



(40) 
(41) 



We let the resonator evolve freely for a time r, and the 
state then becomes 



\$n(t)) = exp(-War)|$Ar(0)), 



(42) 



= -j=\gg)[\0) + exp(-i<p If )\N)], (43) 



where ip?{ is the accumulated phase factor, 



(44) 



Now we study how to transfer the relative phase be- 
tween the two Fock states to the qubit 1. We switch-on 
the interaction between the two qubits and the resonator 
sequentially, for the times t* N and t* N _ x given by 



N 



and t 



N-l 



(45) 



2gWN J " ± 2g l y r N~^l' 

where we set g\ s» gi- The state thus becomes 

|<Mr 2 )> a _L[| 5ff )|o)-exp(-i^)|ee)|Ar-2)], (46) 



where r 2 = r + t* N + t* N _ 1 . We have assumed that the 
time duration r is much greater than the times and 



t 



N-l 



To perform the CNOT gate, we detune the qubits 
from the resonator and set the two qubits at the same 
frequency. We then apply the CNOT gate such that the 
state |ee) will change to \eg). The state then becomes 



l<Mr 2 )) 



1 



_[| 5ff )|0)-exp(- 



^)|e 5 )|iV-2)], (47) 



where t' 2 = r 2 + t c , and t c is the time duration for 
the CNOT gate. We only turn-on the interaction be- 
tween the qubit 2 and the resonator for a time t* N _ 2 = 
ir/2gi\fN — 2, and this gives 

|^(7a))«-L [\ gg )\Q) + i ex p (-i VN )\ee) \N - 3)] , (48) 



where T3 



T'+t 



N-2- 



In the interaction picture, we 

(2), 



repeatedly apply the evolution operator \U 2 U c ) until 
the number of photons in the resonator becomes zero. 
We summarize this procedure in the interaction picture 

as 
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uM 2 \tl)U c f/| 2) (^- 2 )^c^ 2) (^-i)^ (1) (^)l*Jv(T)), 



(49) 



2(JV-2) terms 



where the time t* is 



2g 1 yjn 1 



(50) 



and n is a positive number for n = 1, . . . , N. Qubit 1 now 
completely disentangles from qubit 2 and the resonator. 
The final state becomes 



\9) + exp ^- ihpN - —)\ l e )|l.9>l°>' 

(51) 

where tf = r + td and td is the total time for disentangle- 
ment. We assume that the free-evolution time r is much 
larger than the time td so that we can ignore the relative 
phase accumulated during the time td- Note that the 
relative phase ip^ is encoded on the qubit 1. 

Afterwards, we apply the tt/2 pulse to the qubit 1 and 
then measure the excited state of qubit 1. We can de- 
termine the phase factor from the probability P e of the 
excited state of qubit 1. For simplicity, we set </>i = in 
Eq. ([TT]) . The probabilities of the excited state of qubit 
1 then become 



p(l) p 


4c 


- sin^jv) 


for iV = 


4fc, 




p(2) P 


4c 


+ cos lpn) 


for N = 


4fc- 


1, 


p(3) „ 
J e 


4c 


+ s'mtf N ) 


for N = 


4fc- 


2- 


P (4) „ 


4c 


— COS if n) 


for TV = 


4fc- 


3, 



(52) 



where fe = 1, 



1. Imperfect CNOT gate operations 



the density matrix pf 



r N-l 



Pf 



■ exp 



■[\g)(g\ + \e)(e\+e^p{iip N )\g)(e\ 

{-i<PN)\e){g\]®\g){g\0)(0\, 



(53) 



where tpN — ujNt. Here we have taken the leading or- 
der approximation of the density matrix (see appendix 
B). We notice that qubit 1 cannot be fully disentangled 
from qubit 2 and the resonator by using imperfect CNOT 
gates. We now apply a tt/2 pulse to qubit 1 and then 
measure the excited state of qubit 1. For simplicity, we 
set 0i = in Eq. (fTTj) . The probabilities of the excited 
states of qubit 1 are 



p(l) p 


4c 




- x sin<p N ) 


for N 


= 4fc, 




p(2) p 

e 


4c 


+ e N - 


COS (fix) 


for TV 


= 4fc- 


1, 


p(3) p 


4c 


+ e N - 




for TV 


= 4fc- 


2, 


p(4) p 


4c 


-e N ~ 


COS <^jv) 


for N 


= 4fc- 


3, 



(54) 



where k 



= I,- 



. . , AT/4. The coherence factors, 
e" "smtpjv or e J " _ 1 cos (fiN in Eq. (|54"|) . contain a coeffi- 
cient e^ -1 which is smaller than one. This means that 
the imperfect CNOT gate operations lead to dephasing 
of the qubit 1. 



2. Dissipation in the photon field 

We now take into account the dissipation effect of the 
photon field in Eq. (|3T))) during the free evolution. After 
a time r, the density matrix p of the photon field can be 
written as [35[ 



In realistic situations, the CNOT gate is not perfect 
due to decoherence or experimental constraints. Let us 
briefly examine quantum phase measurements using im- 
perfect CNOT gates. Here we assume that the fidelity 
of this imperfect CNOT gate is very high. The small 
imperfections of this CNOT gate can be characterized 
by a parameter e which is positive and close to one. If 
this CNOT gate would be perfect, then the parameter e 
would be equal to one. A more detailed discussion of the 
effects of non-ideal CNOT gate operations on the disen- 
tanglement process is given in appendix B. We repeat the 
same procedure in Eq. (|49[) using a number of imperfect 
CNOT gate operations. The system can be described by 



P=1 2 



■exp(iip N -TNt)\0)(N\ 



-exp (~iip N ~TNt)\N)(0\ 

V iV!exp(-2fcrr)[l -exp(-2rr)] JV - fe 



k=0 



k\(N-k)\ 



\k)(k\ 



(55) 



We have also assumed that the decoherence of the qubit 
is negligible during the disentanglement process and the 
time td is much smaller than the dissipation timescale 
(rN)^ 1 . The probabilities of the excited state of qubit 1 
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then become 



P e 1] ~ 2( 1+ /") fOT N = 4fc ' 

p e 2) ~ \{ l + h +) fo r N = 4k — l, 

Pi 3 ' » ^(1 + /+) for7V = 4fc-2, 

P e (4) w -(1 + h-) for TV = 4fc - 3, 



(56) 



where k = 1, . . . , iV/4, and /± and /i± are two functions 
given by 

f± = ±e N - 1 exp(-TNr)sm.(p N , (57) 
h± =±e N - 1 exp(-TNT)coscp N . (58) 

3. Kerr nonlinearity 

We note that the superposition state in Eq. JT]) can also 
be used to measure the phase due to the Kerr nonlinearity 
[2 EH HI1 • The Hamiltonian of the nonlinear interaction 
is given by |24| 



TT ^/n^ 2 n 2 

nonlinear — A w u i 



(59) 



where x is the interaction strength. The two coupled 
CPB qubits can form a nonlinear medium of the photon 
field [24| (see also Ref. [HI). The strength \ of the non- 
linear interaction in Eq. (|5^|) can attain 1 GHz or even 
higher 24j]. We only turn-on the nonlinear interaction 
during the free evolution. Thus, the phase factor if>N can 
be rewritten as 



ujNt, 



where 



u = lu + x(N - 1) 



(60) 



(61) 



is an effective frequency due to the Kerr nonlinearity. We 
then apply the same procedure to detect the phase of the 
photon field. 



4- Visibility 

Now we investigate the visibility of the quantum co- 
herence. The visibilities and denote the odd 
and even photon numbers, respectively. From Eqs. (|33[) 
and (f56|) . the visibility V^ 1 ' for an odd number of photons 
can be found as 



jv-i r 



Ti = — < arccos 

UJ 



z 2 + r 2 



exp (-IYi) 



, (62) 
(63) 



where e < 1 is the parameter that characterizes the 
quality of the imperfect CNOT gate and w = cj + 
x(N — 1) is the effective frequency in Eq. (|6Tj) . 

For an even number of photons, we take the form of 
the coherence factor 



CM 



JV-l 



exp (— TNt) sinwr. 



(64) 



The visibility V (2) for even photon numbers can be found 
as (£> > T) 



y(2 



JV-l 



Trr 

exp I -— | +exp 



3?rr\ 



(65) 



The visibility is slightly higher than V^ 2 \ 

For x = an d e = 1, we find that the visibilities 
have no difference between the cases of single-photon 
state l^i) and multi-photon superposition state |\£jv) in 
Eq. (fTJ). The degree of visibility depends on the ratio of 
uj and r. However, the visibilities of the superposition 
of multi-photon states can be greatly enhanced with the 
nonlinear interaction strength %. For e = 1, the degree of 
visibility depends on the ratio Ui/T = (u) + xN)/T from 
Eq. (|6"Tj) such that higher visibilities can be obtained for 
larger number of photons (N 3> 1). This means that 
the superposition of multi-photon states can show a high 
contrast of interference fringes, even in the presence of 
dissipation. 



IV. FACTORIZING INTEGERS USING 
SUPERCONDUCTING CIRCUITS AND GAUSS 
SUMS 

We have shown that the quantum phase of the photon 
states in a superconducting resonator can be measured 
with qubits. These measurement methods are useful for 
the implementation of the Gauss sum (this section) and 
metrology (the following section). 

Now we study a physical realization of the Gauss sum 
algorithm using superconducting circuits. The Gauss 
sum can verify whether a number n is a factor of another 
number N or not. The Gauss sum [H, H3| is defined by 



1 / N 

Cjv(") = ~ cxp ( _ 27r?fc 2 — 



k=0 



(66) 



If n is a factor of N, then |Cjy-(fi)| is equal to one. Oth- 
erwise, the value of |C»v(n)| is less than one. To find 
the factors n of the integer N, the trial factor n scans 

through all numbers from 2 to y/N [U |M SH . Thus, 
in this manner, it can factorize integers. However, this 
Gauss sum algorithm does not provide a speed-up over 
classical computation. 

To save considerable experimental resources and mini- 
mize the effects of decoherence, it is advantageous to use 
as few terms as possible in the Gauss sum in Eq. (|66[) 
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31] 
31] 



Thus, the truncated Gauss sum can be employed 



7, 11,13 



c£(fi) = 



K 



K 
fc=0 



exp 



2mk 2 — 
n 



(67) 



where K is a, positive integer which is smaller than h. We 
only need to sum over K + 1 terms instead of the total 
n terms so that considerable experimental resources can 
be saved. 



A. Single-photon case 

We can apply the same method in determining the 
phase factor of the photon field in the resonator using 
the superposition of the vacuum and single-photon state 
in Eq. ([T]). We follow the same method in Sec. III. A to 
generate the superposition of the vacuum and the single- 
photon state with equal weights, i.e., 



1 



|$ 1 ) = - 7 =| 5 )(|0) + |1)). 



Waiting the following times 



T k = 2irk 2 ^-, 
nu> 



(68) 



(69) 



for the free-evolution, allows the state to accumulate a 
relative phase (—iuTk) between the two Fock states |0) 
and |1). The state then becomes 

|*x (t*)) = -j=\g) [|0) + eM-i"r k )\l)] . (70) 

Using the method described in Sec. III. A, we can trans- 
fer the relative phase to qubit 1. The state of qubit 1 
should then be measured after applying a 7r/2-pulse to 
it. The phase factor cos(u>rk) can be determined from the 
probability of the excited state of qubit 1 in Eq. (1291) . 

By repeating the same procedure K + 1 times and set- 
ting the time duration Tk for k — 0, . . . , K , we take the 
average of the total sum. We readily obtain the real part 
of the truncated Gauss sum in Eq. ([rJTl) as [3l| 



R%{h) 



K 



1 K 

+ 1 ^ 

fc=0 



cos [ 27rfc 2 



N 



(71) 



Indeed, the real part of the truncated Gauss sums has 
been shown to experimentally find the factors of integers 

a 

Now we study how to apply the truncated Gauss sum 
for checking factors. In Fig. [31 the truncated Gauss 
sum \Rf(n)\ for N = 1001 = 7 x 11 x 13 is plotted 
against the trial factors n. The truncated Gauss sums are 
represented by black circles. We have only used K +1 = 5 
terms in the truncated Gauss sum for n by scanning fi 
from 2 to VlOOl w 31. 



0.8 



as 



0.6 



- 

X 



0.4 



0.2 ( h 



10 15 _ 20 25 

n 



30 



FIG. 3. (Color online) Truncated Gauss sum for N = 
1001 = 7 x 11 x 13, plotted versus the trial factors h. Re- 
sults for a damping rate T = (no dissipation) are shown 
with black circles; while T = 6.92 x 10~ 6 to are shown with 
red crosses. We have summed over K + 1 = 5 terms for 
n = 2, . . ., vTOOl w 32. The factors of N are marked by the 
vertical lines located at fi = 7, 11 and 13. All factors are above 
the threshold ~ 0.7071, indicated by a thin horizontal 

dotted line. 



In Ref. [42|, Stefahak et al. gave a threshold to dis- 
criminate factors from non-factors. They found that the 
truncated Gauss sum for non-factors is bounded from 
above by l/y/2 in the limit of large K [42]. As shown in 
Fig. [31 we can see that all factors are above the threshold 
of the Gauss sums. l/y/2 w 0.7071 [13, whereas all non- 
factors are below the threshold, l/y/2 w 0.7071, shown 
by the horizontal dotted line. Therefore, this enables us 
to clearly distinguish the factors from nonfactors. 



1. Effect of dissipation of the photon field 

We now examine the truncated Gauss sum in the pres- 
ence of dissipation of the photon field. We adopt the 
same method as discussed in the previous subsection to 
determine the phase factor. From Eq. ([32]) . the truncated 
Gauss sums can be written as 



R%{h) 



K 



1 K 

k=0 



-Tri, 



2-Kk 1 



,N 



(72) 



The performance of the Gauss sums can decrease due to 
dissipation. The terms with higher k can become vanish- 
ingly small in Eq. (f72"j) . This limits the size of the number 
N to be vertified by the Gauss sum. In Fig. [3l we plot 
the truncated Gauss sum R*~ (n) as a function of the trial 
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numbers h for N = 1001 = 7 x 11 x 13. Here we consider 
the damping rate T/ui = 6.92 x 10~ 6 (we have taken the 
values of w/2tt = 6.57 GHz and T = 2.86 x 10 5 Hz from 
the experiment in Q). We then sum over K+l = 5 terms 
for checking the 31 numbers n — 2, . . . , vlOOl ~ 32. As 
shown in Fig. [3J the truncated Gauss sums slightly de- 
crease due to damping (denoted by the red cross). How- 
ever, we can see that the truncated Gauss sum can still 
be used for distinguishing the factors from nonfactors if 
the relaxation times of resonators are about several /is 
1- 



B. Multi-photon case 
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FIG. 4. (Color online) The truncated Gauss sum with K + 
1 = 5 terms for N = 5005 = 5x7xllxl3is plotted versus the 
trial factors h, for the nonlinear Kerr interaction strength \ — 
u> and damping rate Y — 6.92 x 10~ 6 tt;. The blue diamonds, 
red crosses and black circles represent the different numbers of 
photons: iV = 1, 3 and 5, respectively. The factors of N (here: 
5, 7, 11, 13, 35, 55 and 65) are indicated by the vertical lines. 
All factors are above the threshold l/\/2 « 0.7071, indicated 
by a thin horizontal dotted line. 

Next we consider the multi-photon superposition states 
l^jv) in Eq. ([T]) for the truncated Gauss sum. We first 
produce the superposition of the vacuum and the multi- 
photon state with equal weights, i.e., 



1 



\* N } = -j=\gg)(\0) + \N)). 
We let the system freely evolve for the times 



N 



ffe = 2nk — — 



(73) 



(74) 



where Q — lu + x(N — 1) is the effective frequency in 
Eq. (|61j) if the Kerr nonlinearity is used. Then, a relative 
phase between the two Fock states is acquired. The state 
becomes 

\$ N (f k )) = i= [|0) + exp(-iQh)\N)] . (75) 

We now apply the phase detection method described in 
Sec. III.B for the multi-photon case. After applying a 
7r/2-pulse to qubit 1, we measure the state of qubit 1. 
Then, for an odd number of photons, we can determine 
the phase factor cos(wTfc) from the probabilities of qubit 
1 in Eq. (j5~4")l . The truncated Gauss sum is given by 



R%(n) 



C N-1 



K 



K 
fe=0 



27rfc" 



,N 



(76) 



where e is a parameter originating from imperfect 
CNOT gate operations. In the presence of dissipation of 
the photon field, we can determine the phase factor from 
the probabilities in Eq. (I5rj)) . The truncated Gauss sum 
becomes 



f JV-i A / N\ 

R%(n) = —— £ exp (-TNf k ) cos 27rfc 2 -- . (77) 

k=0 ^ U ' 



C. Kerr nonlinearity 

Now we study the performance of the truncated Gauss 
sum for checking factors using the Kerr nonlinearity. In 
Fig. HI the truncated Gauss sum is plotted versus the 
trial factors n for N = 5005 = 5 x 7 x 11 x 13. The 
truncated Gauss sums with various numbers of photons 
[N = 1 (blue diamonds), N — 3 (red crosses) and N = 5 
(black circles)] are shown in Fig. [4] We can see that the 
truncated Gauss sum for the factors are much closer to 
unity if the larger photon numbers are used. The trun- 
cated Gauss sum using the Kerr nonlinearity can show 
a clearer pattern to discriminate the factors from non- 
factors and a much larger number can be verified [43| . 
This can be easily understood by rewriting the truncated 
Gauss sum [from Eqs. ([71)1 and (|77|)] as 



R%(h) = 



1 



K 



K 



1 ^ 

fc=0 



exp 



27rfc 2 rA> 



cos I 2irk — 

h 



The exponential functions in Eq. (|78|) are closer to one 
when the effective frequency Co = lu + x(iV — 1) becomes 
higher. The larger number of photons leads to a higher 
value of <D. Therefore, the use of Kerr nonlinearities en- 
hances the performance of the truncated Gauss sum for 
finding the factors of an integer, even in the presence of 
dissipation in the photon field. 

We can roughly estimate the size of the factorized num- 
ber N for which the exponential factor in Eq. ([75)1 is 
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around exp (—1). In this case, the fully factorizable num- 
bers are then multiplied by a factor of N and it can attain 
10 2 iV if the nonlinear strength x is about the frequency 
of the resonator, and with the same values for the pa- 
rameters discussed above and K < 10. If the relaxation 
time of the resonator were to be ~ 1 [is, ~ 1 ms and 
~ 1 s, then the largest partially-factorizable numbers 
would be - 10 4 N, ~ l0 7 iV and - 10 10 N, respectively, 
where N is the number of photons in the resonator. 



1. Effect of imperfect CNOT gate operations 
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FIG. 5. (Color online) Truncated Gauss sum with K + 1 = 5 
terms for N = 5005 = 5x7x11x13 plotted versus the 
trial factors h, where the nonlinear Kerr interaction strength 
X = ui, the damping rate F = 6.92 x 10 _6 cj, and the number of 
photons N is 3. Different values of e are shown: e = 1 (black 
circles), e = 0.99 (blue squares) and e = 0.95 (red diamonds). 
The factors of N are indicated by the vertical lines. All factors 
are above the threshold l/\/2 « 0.7071, indicated by a thin 
horizontal dotted line. 

Note that we require more CNOT gate operations to 
disentangle qubit 1 from both qubit 2 and the resonator. 
This requirement is for the multi-photon states involving 
higher number of photons in Eq. ([T]). From Eq. (|77|) . the 
truncated Gauss sum scales with the coefficient e N_1 . 
Thus, the imperfect CNOT gates unavoidably affect the 
performance of the truncated Gauss sum. 

In Fig. [5J we plot the truncated Gauss sum for N = 
5005 = 5x7x11x13 against the trial factors n, using the 
multi-photon superposition state for N = 3 in Eq. ([1}. 
We find that the performance of the truncated Gauss 
sum is very sensitive to small variations of e. As shown 
in Fig. [5J the values of the truncated Gauss sum for 
e = 0.95 (red diamonds) are much lower than the case for 



e = 0.99 (blue squares) and e = 1 (black circles). This 
may limit the use of multi-photon states \^n) for the 
truncated Gauss sum if the photon number N becomes 
large. 



V. PRECISION MEASUREMENT OF THE 
RESONATOR'S FREQUENCY 

We can determine the quantum phase factor of the 
photon states by detecting the state of the qubit. This 
should enable to precisely measure the frequency of the 
resonator ui. The uncertainty |djo;| of the frequency ui of 
the resonator is given by [44 - l4q [ 



\Su 



1 y/P e (l-P e ) 





dP e 











(79) 



where M is the number of measurements. 



A. Single-photon case 

We now consider the detection with the superposi- 
tion of the vacuum and the single-photon state. Using 
Eq. ([29)1. the minimum value of the uncertainty |5o;| m i n , 
for T = 0, is 



\6u\t 



Mr 



(80) 



where r = twk jlio is the time duration of each measure- 
ment and m is an odd integer. 

In the presence of dissipation (r^0), the uncertainty 
is given by 



IM 



1 — exp (— 2r-r) cos 2 (wr) 



Mr 2 exp (-2Tt) sin 2 (wr) 



1/2 



(81) 



We assume that the frequency lj is much greater than 
the damping rate T. The minimum uncertainty |5a;| m i n 
can be found as 



\5u\ 



for lot = rrnr/2 and 



1 



exp(l)r 



T = —■ 



(82) 



(83) 



The precise measurement of the resonator's frequency 
ui can be used to determine the length of the resonator. 
The frequency of the resonator u> is proportional to nc/L, 
where c is the speed of light and L is the length of 
resonator. This enables us to measure the macroscopic 
quantity L precisely and its accuracy is up to 



|*L| = 



L\5ui\ 



UJ 



(84) 
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B. Multi-photon case 



ACKNOWLEDGMENTS 



Next we study the degree of accuracy in the phase 
measurement with the superposition of the multi-photon 
state |*jv) in Eq. {TJ. Without dissipation (T = 0), 
the minimum value of the uncertainty |<5w| m i n of the fre- 
quency lu is 



IM 



1 



M Nt 



(85) 



This accuracy scales as e 1 /N with the superposition 

state | \E r jv) which gives a much better improvement than 
using the single-photon state j^i) if e is very close to 
one. However, in the presence of dissipation (r^O), the 
minimum uncertainty of |5w| m i n becomes 

1 



exp(l)e 



l-Ni 



(86) 



It is no different to the single-photon case even if e = 1. 

However, the nonlinear interaction strength \ can be 
measured precisely. The minimum uncertainty |<$x|min is 
given by 



\SxU 



1 exp(l)e 1 ~ A 'r 



N 



M 



(87) 



This uncertainty scales with e 1_Ar /N under decoherence. 
This enables to detect the strength of the Kerr nonlin- 
ear ity with very high accuracy if e 1 ^ 1 ^ /N is very small. 



VI. SUMMARY 



We have presented a method to measure the quantum 
coherence of the superposition of two number states in 
a superconducting resonator with two Josephson phase 
qubits. We have also studied the visibility of the quantum 
coherence in a dissipative environment. We found that 
the visibility of the photon can be enhanced if nonlinear 
interactions are used. This may be useful to probe the 
quantum coherence of multi-photon superposition states. 
We showed that the phase measurement scheme can be 
applied to factorizing integers and parameter estimation. 

The detection of the superposition of the vacuum and 
the single-photon state can be realized with current tech- 
nology [3]. But the measurement of multi-photon super- 
position states involves a number of CNOT gate opera- 
tions to disentangle the qubit from the resonator. The 
quality of measurements is degraded due to the imperfect 
operations of the CNOT gate. Also, significant resources 
can be consumed when measuring the superposition of 
multi-photon states because the number of gate opera- 
tions is proportional to the number of photons N being 
detected. Here we proposed a method which can be used 
to detect superposition states with a few photons. 
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Appendix A: Qubit-resonator coupling in the 
far-detuning regime 

The Hamiltonian ffO) = H les + H% + H%_ Ies 

— ua^a + ""rp^jz + 9j( a<T j+ + °j- a *) (Al) 

can be exactly solved [23| for j = 1,2. We now consider 
the subspace spanned by the basis |<7}|n) and \e)\n — 1). 
The two eigenvalues can be solved as 



and the corresponding eigenvectors are 



= sm^|s)|n) + cos^|e)|n-l), (A3) 
I-),- = coa pj\g)\n) - single)!?-! - 1), (A4) 



where 



sin/3j = — 

COS/Jj = 



A, + 5j 



v /(A J +^)2 + 4 5 |n' 
2g jy /n 



A, = 



(A5) 
(A6) 

(A7) 
(A8) 



If the detuning Aj = ojQj — uj is much larger than the 
interaction strength gj, the two eigenvectors becomes 



Ie)k-1) 



(A9) 



in Eqs. (|A3[) and (|A4[) . respectively. Therefore, the in- 
teraction between the qubit and the photon field can be 
effectively turned-off by far-detuning the qubit from the 
resonator. 



Appendix B: Analysis of imperfect CNOT gate 
operations on the disentanglement process 

In this appendix, we study the effects of imperfect 
CNOT gate operations on the disentanglement process. 
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We consider the CNOT gate operations to be described 
by a process which suffers from dissipation, decoherence 
and imperfections |47l| . In general, this operation can be 
represented by £ [47] such that 



Pout = £(Pin), 



(Bl) 



where p m and p Q ut are the density matrices of the in- 
put and output states, respectively. The operation £ 
is a convex- linear map and also a positive map [48| . 
Here we assume the operation £ is trace-preserving, i.e., 
tv[£{ Pm )] = l. 

It is convenient to write the states as 



|0) = \gg), |1> = \ge), 
|2) = |ee), |3) = \eg). 

An ideal CNOT gate is defined as 

|omo), |i)->|i>, 

|2>->|3), |3>->|2). 



(B2) 



(B3) 



Now we consider the non-ideal CNOT gate operation but 
with a high fidelity as: 



<0|£(|0>(0|)|0> 
(3|£(|2)(2|)|3) 
(0|£(|0)(1|)|1) 
(0|f(|0)(3|)|2) 
(l|f(|l>(3|)|2) 



(l|£(|l)<l|)|l)*e, 
(2|£(|3)<3|)|2> « e, 
(0|£(|0)(2|)|3)« £ , 
(l|£(|l)<2|)|3)*e, 



(B4) 
(B5) 
(B6) 
(B7) 
(B8) 



and £(\i)(j\)* = £(\j)(i\), where e is a positive num- 
ber close to one. Otherwise, the remaining terms 

(i'\£(\i)(j\)\j') are equal to small parameters e\- , which 
are much smaller than e, where i, i' , j and j' — 0, 1, 2, 3. 



We now follow the same disentanglement procedure as 
summarized in Eq. (|4T)|) , by using the imperfect CNOT 
gates. The density matrix of the total system becomes 



Pf = — [bXfll + |e)(e| + exp{i<p N )\g)(e\ 

+ exp(-i<PN) |e> (g\]®\g) (g\0) (0| + G{< 



.N-2i'j'\ 



(B9) 



where tpN — ujNt. The leading order approximation is 
the density matrix containing the first four terms with 
the coefficients e w_1 . The remaining terms involving a 
large number of entangled states of qubit 1, qubit 2 and 
the resonator are of the order of e . Using imper- 

fect CNOT gates, qubit 1 cannot be completely disen- 
tangled from qubit 2 and the resonator in Eq. (|B9I) . 
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